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Thermodynamic  Extrapolation  of  Rocket  Engine 
Performance  Parameters 
Ge  Minglong 

Abstract 

In  this  paper,  an  isentropic  reference  line  and  two  partial 
derivatives  with  fixed  enthalpy  were  used  to  establish  the  ex¬ 
trapolation  formulas  for  rocket  engine  performance  parameters. 
When  the  initial  enthalpy  value  of  the  propellant,  combustion 
chamber  pressure,  and  nozzle  exit  pressure  began  to  vary,  these 
formulas  could  be  used  to  carry  out  extrapolating  calculations 
for  specific  thrust,  characteristic  velocity,  nozzle  area  ratio, 
combustion  chamber  thermodynamic  parameters,  and  nozzle  exit 
thermodynamic  parameters.  The  formulas  in  this  paper  are  simple 
and  easy  to  use.  The  accuracy  is  high. 

I.  Introduction 

The  thermodynamic  calculation  of  the  theoretical  performance 

characteristics  of  a  rocket  engine  is  a  classical  case.  It  has 

been  introduced  in  Reference  1  as  well  as  in  other  books  and 

journals  concerning  rocket  engines.  This  complicated  calculation 

is  very  useful  in  the  design  of  a  rocket  engine.  However,  it 

consumes  a  considerable  amount  of  computer  time.  Furthermore, 

it  is  also  inconvenient  to  organize  the  great  deat  of  data  into 

volumes.  Therefore,  the  extrapolating  calculation  problem  of 

the  theoretical  performance  characteristics  of  a  rocket  engine 

had  been  brought  out.  Some  people  gave  the  formulas  in  thermo- 
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dynamic  extrapolating  calculation. 

In  recent  years,  along  with  the  availability  of  pocket 
electronic  calculators,  it  is  even  more  convenient  to  use  an 
extrapolation  method  to  calculate  the  theoretical  characteris¬ 
tics  of  a  rocket  engine.  For  this  reason,  on  the  basis  of 
presently  available  extrapolation  formulas,  more  accurate  extra¬ 
polation  formulas  were  derived  in  this  paper. 

II.  Isentropic  Reference  Line 

For  a  certain  propellant  composition  at  a  fixed  mixture 
ratio  0/F,  because  its  elemental  composition  is  fixed,  therefore, 


1 


the  ideal  thermodynamic  process  of  the  engine  -  the  isentropic 

limiting  equilibrated  expansion  process  -  can  be  expressed  as 

various  isentropic  lines  such  as  c  t  e  e",  cte,...etc.,  on  the 

o  o  o 

same  isentropic  diagram  as  shown  in  Figure  1  regardless  of  the 

same  isentropic  value  (Hcq  or  Hc)  of  the  propellant  as  affected 

by  factors  such  as  the  environmental  temperature.  It  is  also 

irrespective  to  the  combustion  chamber  pressure  (pco»  pc' '  or  ' 

and  the  nozzle  exit  pressure  (P  ,  P  1 ,  P  "  or  P  ) . 

r  co'  e  '  e  e 


Figure  1.  Schematic  Diagram  of  Enthalpy  and  Entropy 

1.  enthalpy  H 

2.  entropy  S 

3.  isentropic  reference  line 

For  a  given  propellant  initial  enthalpy  value  (Hcq)  an^ 
combustion  chamber  pressure  (PCQ) •  it  is  possible  to  determine 
the  thermodynamic  parameters  for  a  combustion  chamber  condition 
point  cQ  and  throat  condition  point  tQ  through  an  accurate 
thermodynamic  calculation.  Then,  several  nozzle  exit  pressure 
levels  were  chosen  to  calculate  the  thermodynamic  parameters  of 
these  nozzle  exit  condition  points  (eQ,  e" ,...). 
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Table  1.  Parameters  on  the  Isentropic  Reference  Line 


1.  (hydrogen  gas  and  oxygen  gas  propellant,  mixing  ratio 
0/F=7.9365,  initial  enthalpy  value  H  =0  Kcal/kg) 

2.  parameter  c 

3.  combustion  chamber 

4 .  throat 

5.  nozzle  exit 

6.  P,  kg/cm^  (absolute) 

7.  H,  Kcal/kg 

8.  M,  kg/kg  mole 

9.  Cp  Kcal/kg  K 

10.  *data  obtained  through  extrapolation. 
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On  the  isentropic  line  c0t0e0e"»  there  are  several  points  with 

known  values  of  sych  parameters  as  pressure  ratio  pc/p,  pressure 

p,  temperature  T,  enthalpy  H,  molecular  weight  M,  and  constant 

pressure  specific  heat  cp,  as  well  as  partial  derivatives  such 

as  /  alnM\  and  /  alnM 
yainp  [  ainT 

two  partial  derivatives  have  been  shown  in  References  1  and  4). 
Therefore,  this  isentropic  line  can  be  used  as  a  reference  line 
in  the  extrapolation.  The  parameters  of  a  typical  isentropic 
reference  line  are  listed  in  Table  1 .  The  parameters  in  the 
upper  half  of  the  table  came  from  Reference  5. 

On  an  isentropic  line,  the  thermodynamic  parameters  of  a 
point  other  than  a  parametric  point  can  be  calculated  by  using 
the  following  extrapolation  formula  in  approximation: 

F  —  £Hiy  | 

(AF),- P(Alnf),  J  (1) 


)(the  calculation  formulas  and  methods  of 
P 


where  F  repreents  parameters  such  as  T,  H,  M...etc.  Q  is  the 
extrapolation  coefficient,  which  is  determination  by  the  para¬ 
meters  of  the  two  known  neighboring  points. 

With  respect  to  the  isentropic  reference  line  in  Table  1, 
the  results  of  the  calculation  of  the  extrapolation  coefficients 
such  as  Qm,  Q„ ,  and  Q„  are  listed  in  the  lower  half  of  the  table. 

I  M  n 

III.  Isenthalpic  Partial  Derivatives  and  Isenthalpic  Correla¬ 
tions 
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Referring  to  the  relevant  formulas  in  References  4  and  5, 
the  two  isenthalpic  partial  derivatives  used  in  this  paper  can 
be  expressed  by  using  the  following  formulas: 


(2) 


(3) 


where  R  is  the  commonly  used  gas  constant  (1.98726  Kcal/kg  mole  K) . 

With  respect  to  the  isentropic  reference  line  in  Table  1 , 
the  values  of  D^, ,  ,  QDT  and  QDM  calculated  from  the  formulas 

given  above  are  listed  in  the  lower  half  of  the  table. 

The  correlations  of  temperature  and  molecular  weight  on  an 
isentropic  line  are 


(Ain T)«  “  (4) 


(AlnM)*  —  D^Alnj)*  (5) 


Substituting  the  expression  of  the  first  law  of  thermo¬ 
dynamics  and  the  equation  of  state  into  the  differential  equa¬ 
tion  of  entropy,  one  gets 


di 


dH  _R  df 
T  U  t 


(A) 
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For  an  isenthalpic  line  (dH=0) ,  the  formula  gives 


(A/)*  — 


(6) 


IV.  Establishment  of  Extrapolation  Formulas 

1 .  Combustion  Chamber  Thermodynamic  Parameters  and  Formula 
for  Characteristic  Velocity. 

When  the  initial  enthalpy  of  the  propellant  Hc^  Hcq  and  the 
combustion  chamber  pressure  pc/  p  ,  based  on  AH  =  H  -H  and 
the  known  data  at  point  cq,  it  is  possible  to  calculate  the 
parameters  of  point  c'  in  Figure  1  using  equation  (1).  Then, 
according  to  the  parameters  at  point  c'  and  p=Pc,  the  parameters 
at  point  c  can  be  calculated  according  to  equations  (4)  and  (5) . 
Finally,  the  combustion  chamber  temperature  and  moelcular  formulas 
can  be  obtained: 


T'~ 


V,- 


( 


(r„  + 


(Dtco  +  I 

fb-&-. 

Qm  ‘  * 

'  P" 

Qh  /J 

(  Dyco  +  QjU&Hi' 

)(b 

1 

Qm  1  L 

'  Qh  J 

'  P" 

Qh>  J 

(7)  /326 

(8) 


Since  the  characteristic  velocity  c*  is  proprotional  to 
v/Tc/Mc ,  the  extrapolation  formula  for  characteristic  velocity 
can  be  written  as 


*•  -  ttJTMJM.T'. 


(9) 


2.  Formulas  of  Thermodynamic  Parameters  at  the  Nozzle 
Exit  (p=Pc) 
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Applying  equation  (6)  between  points  c  and  c' ,  and  e  and  e' 


*  H,  f. 


(A) 


Substitute  equation  (1)  into  this  formula,  one  gets 


In  IjL  -  in_f£.  + 

Pro  Pro 


u 


co 


Qu±H, 

Qh 


(10) 


where  Meo  is  the  approximate  value  of  Me' .  It  is  the  known 
molecular  weight  of  p  =  peo  on  the  isentropic  reference  line 
close  to  pe' .  In  order  to  select  p  ^  so  that  M  _  can  be  deter- 
mined,  the  initial  value  of  pe'  can  be  calculated  according  to 
the  following  formula: 


po-m 


(11) 


This  formula  was  obtained  as  an  approximation  of  equation  (10) . 
The  Me"  in  the  equation  is  the  known  molecular  weight  of  the 
lowest  point  on  the  isentropic  reference  point. 

The  values  near  peQ  were  adopted  as  the  extrapolation  co¬ 
efficients  Qt,  Qdt,  Qm,  Qdm,  Q  ,  .in  .equations  (7),  (B),  (10)  and 
(11) .  When  AHc<0,  the  average  value  of  cQ-t  was  used.  When 

AH  >0 ,  the  values  for  p>p  were  used, 
c  rco 

Based  on  the  known  data  at  point  e  and  the  value  of 
I  P  ° 

In  _e  calculated  from  equation  (10)  ,  the  parameters  at  point  e' 
Peo 

can  be  calculated  using  equation  (1).  Then,  based  on  the  par¬ 
ameters  at  point  e'  and  P=P  ,  the  parameters  at  point  e  can  be 

V 

calculated  from  equations  (4)  and  (5).  Thus,  the  formulas  for 
nozzle  exit  temperature  and  molecular  weight  were  obtained. 
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T,  -  (t..  +  Qrin  «p  [(  Dr,.  +  £>Drln-^)(ln  -  In  -J^)] 

M,  —  (a/,.  ■+*  i>«in  ctp  [Yd*,.  +  Po^ln-^-Vln  ~^£-  In 

'  '  *■  '  Pf  '  '  £»•  P#o  J 


3.  Formulas  for  parameters  such  as  the  specific  Thrust 

Substituting  H  =  H  .  =r  H  +  AH  and  H  =  H  .  =  H  +  Q„ 

3  c  c*  co  c  e  e  *  eo  H 

|n  Pe ’  into  the  specific  thrust  formula  at  the  designed  altitude 
P 

eo 

I  =  294.98  /(Hc-He) /1000,  we  get 


/-  9.328  JHt.  +  £iH.-H..~QH\n-^L  (14) 

*  P •• 

The  extrapolating  coefficients  ,  QDT ,  QM,  QDfJ,  and  Qh  in  equa¬ 
tions  (12)  through  (14)  adopted  their  values  at  near  P  .  When 

:  p  e° 

In  e *  <0,  the  values  less  than  those  of  PeQ  were  used.  When 
Peo 

In  P  .  >0,  the  values  greater  than  those  of  P  were  used, 
e '  3  eo 


The  nozzle  exit  area  ratio  and  the  nozzle  area  ratio  are  /327 
calculated  from  the  following  formulas: 


f,  -  *6,50T,/Ir.M. 


(15) 


e  -  UP,/'* 


(16) 
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V.  Accuracy  of  Extrapolation  Formulas 

On  the  basis  of  the  formulas  given  above,  it  is  possible  to 
carry  out  extrapolation  calculations.  For  a  hydrogen  oxygen 
propellant  with  a  mixing  ratio  0/F=7.9365,  some  of  the  results 
of  extrapolation,  such  as  nozzle  area  ratio,  characteristic 
velocity,  combustion  chamber  temperature,  and  specific  thrust  at 
designed  altitude,  are  listed  in  Tables  2,  3,4  and  Figure  2, 
respectively.  For  comparison,  the  results Q f  t^e  correspond¬ 
ing  accurate  calculation  and  the  extrapolated  calculation  in 
Reference  5  are  also  listed  in  these  figures  and  tables. 


f.h 

j.  *««*,**  H«- -190.4  *+/£ff 

Z.*l*ft*,»*  H.- 629.1  **/£Jf? 

j.  Una 

I,***)*** 

a* s  ft*** 

( ■  »«** 

7**ft**« 

iesftft** 

10 

2.4(8 

2.47( 

2.289 

2.4(6 

2.459 

2.337 

40.83 

7.151 

7.188 

7.130 

7.169 

7.196 

7.279 

«.0S 

10.75 

10.80 

10.77 

10.81 

10.77 

11.00 

too 

14.(9 

14.74 

14.70 

14.81 

14.70 

15.01 

300 

36. n 

3(.40 

35.70 

37.02 

37.05 

36.45 

400 

46.04 

4(.19 

45.04 

47.22 

47.29 

45.98 

(00 

(4.44 

(4.58 

(2.49 

((.(6 

66.55 

(3.80 

•00 

•1.77 

•1.90 

71.83 

•5.21 

•4.97 

•0.49 

1000 

98.31 

98.41 

94.40 

103.2 

103.6 

96.38 

Table  2.  Comparison  of  the  Nozzle  Area  Ratio  Calculated  by 
Extrapolation  with  Its  Accurate  Value 

(hydrogen  oxygen  propellant,  mixing  ratio  0/F=7.9365, 
combustion  chamber  pressure  p  =6. 805Kg/cm2 ,  absolute) 

C 

1.  liquid  hydrogen,  liquid  oxygen,  enthalpy  H  = 
-190.6  Kcal/kg 

2.  gas  hydrogen,  ozone,  enthalpy  Hc  =  629.1  Kcal/kg 

3.  accurate  value 

4.  extrapolated  value  in  this  work 

5.  extrapolated  value  in  Reference  5 

6.  accurate  value 

7.  extrapolated  value  in  this  work 

8.  extrapolated  value  in  Reference  5. 
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l .  tt* «£*  f. 

i  h.  -  o 

i .  *+/ 

**««,&*  H,  - 

3-  -190.6  *-*r/£Fr 

n***»**».- 629.1 
4-.  *+/&ff 

5\ 

HP* 

<•*  *  « 

1  <• 

«*«* 

it  * 

13 

*  *  ft  * 

"*X 

\«5 

68.05 

2197 

2197 

2197 

2156 

2156 

2157 

2324 

2326 

2329 

40.83 

2184 

2182 

2184 

2144 

2142 

2144 

2308 

2310 

2316 

6.805 

2136 

2132 

2139 

2099 

2094 

2099 

2251 

2254 

2271 

n.  tt:  ********«****/* 


Table  3 .  Comparison  of  the  Extrapolated  Characteristic  Velocity 
with  Its  Accurate  Value  (hydrogen-oxygen  propellant, 
mixing  ratio  0/F  =  7.9365) 

2 

1 .  combustion  chamber  pressure  p  ,  Kg/cm  ,  absolute 

2.  gas  hydrogen,  gas  oxygen,  entfialpy  H  =  0  Kcal/kg) 

3.  liquid  hydrogen,  liquid  oxygen,  enthSlpy  H  = 

-190.6  Kcal/kg  c 

4.  gas  hydrogen,  ozone,  enthalpy  H  =  629.1  Kcal/kg 

5.  accurate  value  c 

6.  extrapolated  value 

7.  this  work 

8 .  Reference  5 

9 .  accurate  value 

10 .  extrapolated  value 

11.  this  work 

12.  Reference  5 

13.  accurate  value 

14.  extrapolated  value 

15.  this  work 

16.  Reference  5 

17.  Note:  the  unit  for  the  numerical  value  of  charac¬ 
teristic  velocity  in  the  table  is  m/sec. 
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ttttXSft  t.te 

1.  ff /**',*# 

1.02 

2.04 

4.08 

10.21 

20.42 

40.83 

61.25 

01.56 

2,  MAX 

3039 

3127 

3217 

3341 

3437 

3534 

3591 

3632 

J . 

3053 

3138 

3225 

3344 

3437 

3533 

3590 

3631 

4.  tt:  K 

Table  4.  Comparison  of  Extrapolated  Combustion  Chamber  Tempera¬ 
ture  with  Its  Accurate  Value 

(liquid  hydrogen  -  liquid  oxygen  propellant,  mixing 
ratio  0/F  =  7.9365) 

2 

1.  combustion  chamber  pressure,  pc  Kg/cm  absolute 

2.  accurate  value 

3.  extrapolated  value  in  this  work 

4.  Note:  the  unit  of  combustion  chamber  pressure 
in  the  table  is  K. 

From  Tables  2  through  4  and  Figure  2  one  can  see  that  the 
accuracy  of  the  extrapolation  formulas  in  this  paper  is  rela¬ 
tively  high.  They  are  suited  for  the  extrapolation  of  the  actual 
parameters  of  a  rocket  engine . 


Figure  2.  Comparison  of  Extrapolated  Specific  Thrust  at  Designed 
Altitude  and  Its  Accurate  Value 

(hydrogen-oxygen  propellant,  mixing  ratio  0/ F=7.9365, 
combustion  chamber  pressure  pc  =  6.805  Kg/cm2, 
absolute) 

1.  specific  thrust  at  designed  altitude,  I,  sec. 

2.  gas  hydrogen-ozone  (enthalpy  H  —  629.1  Kcal/kg) 

3.  liquid  hydrogen-liquid  oxygen  Xenthalpy  H  = 

-190.6  Kcal/kg) 

4 .  accurate  value 

5.  extrapolated  value  in  this  work 

6.  extrapolated  value  in  Reference  (5) 

7.  pressure  ratio  Pc/P 
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A  Thermodynamic  Calculation  Method  Of  The  Symbolic 
Formula  For  A  Rocket  Engine 
Fang  Zhaokui 

In  the  past,  a  chemical  formula  thermodynamic  calculation 
was  widely  used  in  our  country.  Its  special  feature  is  that: 
the  chemical  reaction  must  be  correlated  to  the  specific  chemical 
elements  of  the  propellant.  This  makes  it  difficult  to  design  a 
generalized  program.  The  symbolic  thermodynamic  calculation 
presented  in  this  paper  eliminated  the  aforementioned  difficulty. 

I.  Symbolic  Elements  and  Molecular  Formula 


For  a  given  pair  of  positive  integers  i  and  j ,  any  propel¬ 
lant  with  j  elements  and  producing  i  combustion  products  is  de¬ 
fined  by  a  series  of  symbolic  elements  A^  and  the  corresponding 
symbolic  element  atomic  number  a ^ .  The  sequential  arrangement  of 
the  symbolic  elements  A ^ : 


n  'A***  (i  —  i,  2  »•••«) 


would  represent  the  moelcular  formula  of  the  ith  combustion 
product  after  the  propellant  is  combusted.  The  symbol  ^  denotes 
the  sequential  arrangement  as  shown  in  equation  (1) . 

For  example,  in  the  case  of  a  propellant  with  four  elements 
C,  H,  0,  and  N,  the  symbolic  elements  A^,  b-2'  an<*  A4  are  use<* 
to  represent  these  four  elements.  Let  the  i^h  combustion  product 
be  CO.  Obviously,  according  to  the  definition  its  symbolic  ele¬ 
mental  molecular  formula  is: 

CO-  IT  -V*  (2) 


The  sequence  of  i  and  j  and  its  correlation  to  the  products  and 
elements  is  arbitrary.  However,  once  the  correlation  is  estab- 
lsihed,  equation  (2)  uniquely  defines  CO. 


II.  Chemical  Equilibrium  Equation  of  Symbolic  Element  Combus¬ 
tion  Product 

Definition:  For  any  molecular  combustion  product  jj  a^, 
the  dissociation  reaction  correlated  to  the  symbolic  element  A^ 

i  i 

IT  '*/•*—*“  2  O'  —  2* •  •»)  (3) 


is  made  to  be  the  chemical  equilibrium  reaction.  Due  to  the  fact 
that  symbolic  elements  are  given  a  part  of  the  definition  of  a 
number,  therefore,  they  obey  the  conventional  multiplication 


a  .  .  X  A 
Di 


j  “  la  A 


aii=0 


From  the  chemical  equilibrium  reaction  theory,  the  chemical 
equilibrium  equation  corresponding  to  (3)  is: 


'n-w /("*'  x  fX  *  •  ■ •  ■ ’  * C) “ 0 


where  P  A  a  and  Pft  represent  the  partial  pressure  of  the 

symbolic  element  molecular  product  jfjAjaji  and  atomic  product  Aj 
in  the  combustion  process,  respectively.  f^(T)  is  the  equilibria 
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constant  of  the  dissociation  reaction  at  temperature  T.  If  a  /3 
certain  element  does  not  exist  in  this  molecule  at  all  (for  ex¬ 
ample,  CO  molecule  do  not  contain  hydrogen  at  all),  P£=l  due  to 
the  conventional  mathematical  law.  Furthermore,  equation  (5)  is 
satisfied.  Taking  the  logarithm  of  equation  (5)  on  both  sides, 
we  get: 


Using  equation  (6) ,  it  is  possible  to  convert  the  moelcular 
partial  pressure  into  atomic  partial  pressure  to  eliminate  the 
molecular  variables.  Consequently,  the  number  of  equations  is 
greatly  reduced. 


III.  Symbolic  Element  Propellant  Matrix  and  Its  Normalization 
Treatment 

Definition:  Assuming  that  there  is  a  complicated  propel¬ 
lant  comprised  of  n  types  of  pure  chemical  substances ,  the 
molecular  formula,  enthalpy,  and  weight  percent  of  the  mixture 

of  each  pure  chemical  substance  are  A  A. a.  ,  H  ,  and  X  (n=l,2.. 

11  3  jn  n  n 

n) ,  respectively.  Then,  the  characteristics  of  this  propellant 
can  be  determined  by  the  matrix 


A  (»,  /  +  2)  — 


(•u  ' 


•  •n  W,  x, 

*  x/i  W,  x, 

H.  x. 


(7) 


(7)  is  called  the  propellant  matrix.  In  the  performance  calcu¬ 
lation,  it  is  only  necessary  to  consider  the  theoretical  charac¬ 
teristics  possessed  by  a  unit  mass.  For  this  reason: 

Definition:  Choose  the  unit  mole  propellant  symbol  ele¬ 
mental  molecular  formula  whose  molecular  weight  is  1  as  the 
standardized  equivalent  formula  of  the  propellant.  The  process 
of  converting  from  the  initial  propellant  into  the  symbolic 
element  propellant  equivalent  formula  is  called  the  normaliza¬ 
tion  process  of  the  symbolic  element  propellant. 

Theory:  Assuming  that  there  is  a  complex  propellant  (7) , 

Uj  is  the  atomic  weight  of  the  symbolic  element  Aj. 
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The  matrix 


B  (»)- 


(8) 


Then,  the  normalization  of  the  propellant  (7)  is  the  realization 
of  the  following  matrix  operation. 


*(j  +  0  -  /*'(«»/  4-  1)  X  «(«) 


(9) 


where  A* (n, j+1)  is  the  transformation  matrix  of  the  matrix  A 
(n,j+l)  which  was  obtained  by  eliminating  the  last  column  in 
matrix  (7).  (9)  is  a  column  matrix  after  normalization.  The 

subscript  "0"  of  a  column  matrix  SjQ  indicates  that  it  does  not 
belong  to  the  "i"  series  of  the  combustion  chamber.  From  equa¬ 
tion  (9),  one  can  derive 


(10) 


This  is  a  natural  result  of  the  normalization  treatment.  For  a 
dual  component  propellant,  if  the  incendiary  agents  correspond¬ 
ing  to  the  forms  of  (7)  and  (8)  are  A(n,j+2)  and  B(n) ,  the 
oxidants  are  A(m,j+2)  and  B(m),  and  the  mixing  ratio  is  K,  then 
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the  normalization  process  of  this  complex  dual  component  propel¬ 
lant  is:  /332 


AQ+  o-  +  1)  X  +  1)  X  #(■)]/(  1  •+•  K)  (11) 

The  final  matrix  obtained  is  a  column  matrix.  The  proof  for  (9) , 
(10) ,  and  (11)  is  omitted. 

IV.  Basic  Equation  And  Program  Design  In  Symbolic  Formula 
Thermodynamic  Calculation 

Let  us  assume  a  propellant  is  defined  by  positive  integers 
i  and  j  and  the  normalization  matrix  A(j+1).  Then,  the  standard 
form  of  the  combustion  reaction  is: 

U  ri  A** - *  ±  Pi  IT  -V*  (12) 

The  formula  above  shows  that  M  moles  of  propellant  produces  com¬ 
bustion  products  yi  Ajaji  at  Partial  pressure  when  it  combusts. 
However,  when  a  condensed  phase  exists,  P^  should  be  considered 
as  the  mole  number  of  the  condensed  phase  product.  Now,  the 
four  sets  of  basic  equations  in  symbolic  formula  thermodynamic 
calculation  are  re-written  as  the  following: 

1.  Mass  Conservation  Equation:  J-, 

1“  1.  2, - >  (13) 

2.  Energy  Conservation  Equation: 

Hn  =  y  (H®).P./M  (isenthalpic-combustion  chamber  state) 

U  ^  Til 

(14) 

Sq  =  £  t (S$) P^/M  (isentropic-combustion  chamber  state) 
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3.  Chemical  Equilibrium  Equation: 


(Fj),/(RD+lB«i-2^I(fr),/(«T)+lnt.,]-0  .-1,2  (15) 

where  ( H.0 )  ^ ,  (S®)^,  and  (F^)i  in  equations  (14)  and  (15)  are  the 
enthalpy,  entropy,  and  free  energy  of  the  i^*1  product  at  temper¬ 
ature  T,  respectively. 

_  R  nPi  (for  gas  phase)  _  P.  (for  gas  phase) 

fci  n  (for  condensed  phase)  wi  ,  ,  .  . 

0  c  1  (for  condensed  phase) 

R  is  the  gas  constant  of  a  mole  of  generalized  gas.  The  formula 
for  the  condensed  phase  is  omitted. 

4.  Dalton's  Law 

PQ  =  £  Pi  (for  gas  phase  only)  (16) 

Solving  the  aforementioned  equation  with  known  A(n,j+2), 
A(m,j+2),  K,  PQ ,  Hq ,  and  becomes  the  basic  problem  in  thermo¬ 
dynamic  calculation. 

In  designing  the  program,  it  is  necessary  to  linearize  the 
equations.  Furthermore,  (6)  is  used  to  reduce  the  number  of 
equations.  Partial  derivative  approximation  and  convergence 
factor  are  used  to  resolve  the  arbitrariness  of  the  initial  value 
in  a  large  range.  We  used  the  BCY  language  to  compile  the  gen¬ 
eralized  program  for  thermodynamic  calculation.  Several  dozens 
of  propellants  were  calculated  on  a  109-B  computer  with  success. 

Due  to  the  limitation  of  pages,  it  is  not  described  in  detail. 
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